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o : 

We perform dynamical calculations of direct muon-transfer rates from thermalized muonic hy- 
drogen isotopes to bare nuclei C 6+ and 8+ . For these three-body charge-transfer reactions with 
Coulomb interaction in the final state we use two-component integro-differential Faddeev-Hahn-type 
^""^ ' equations in configuration space using close-coupling approximation. To take into account the high 

polarizability of the muonic hydrogen due to the large charge of the incident nuclei, a polarization 
potential is included in the elastic channel. The large final-state Coulomb interaction is explicitly 
included in the transfer channel. The transfer rates so calculated are in good agreement with recent 
experiments. 

B\ 

2 ■ I- INTRODUCTION 

■ Mainly motivated by the possibility of muon-catalyzed fusion of the hydrogen isotopes, theoretical and experimental 
investigations of exotic atomic and molecular systems involving negative muon (/i~) and reactions in such systems 
continue to be active fields of current research pF . Particular attention is devoted to the study of the muon-transfer 
reaction from muonic hydrogen to other elements X z+ with large positive charge Z, because such reactions may have 
| large cross sections and rates. If the hydrogen is contaminated by even a small amount of these heavier elements, 
, this may strongly influence the process of muon-catalyzed fusion by hydrogen isotopes. Consequently, there has 

been considerable recent experimental interest in the study of the muon-transfer reaction in the collision of the 
' muonic hydrogen by heavier nuclei, e.g. carbon (C 6+ ), oxygen (0 8+ ), neon (Ne 10+ ), argon (Ar 18+ ) etc. On 

the theoretical side, these three-body charge transfer reactions involving a heavy transferred particle like muon and a 
qq . strong Coulomb interaction in the final state involving nuclei, such as carbon and oxygen, continue to be challanging 
^C) ' problems calling for new investigations. 

Experimental study of muon-transfer rates from muonic hydrogen H M to heavier nuclei, such as Ar 18+ , Kr 36+ and 
CNJ i Xe 54+ , have revealed a smooth Z dependence M. In these cases the transfer rate Xz increases linearly with Z. 
Theoretical analyses are also in agreement with this conclusion However, it has been found in recent experiment 
that the predicted monotonic Z dependence of the muon-transfer rate is not valid for all Z la|^]. For these transfer 
rates, pronounced fluctuations have been observed for elements up to argon contrary to the smooth Z dependence. 
""^5 ' The experimental muon transfer rates for reactions like 

o : 

•53- (fL J ) ls +X z +-,(X M )( z - 1 )++H+ (1.1) 

>>'... n 

depend in a complicated manner on the charge Z |9l . Here H stands for the hydrogen isotope p (proton) or d (deuteron) 

Q-r and X z+ stands for the target nuclei. 

Another phenomenon which has not yet found a satisfact ory theoretical explanation is the measured isotope effect, 



e.g. the trend of the direct muon-transfer rate of reaction (1.1) from p/j, and d\x to X z+ , where X z+ represents C 6+ , 



^ 8+ || , Ne 10+ @, Ar 18 + [§, and Xe 54+ @. In cases of 0*+, Ar 18 + and Xe 54 + the direct-transfer rate decreases 
with increasing mass of the nydrogen isotope. Theoretical analyses [j3||TTl also support this trend. The experimental 
- ■ results for Ne 10 ||^] and sulphur dioxide differ considerably from the theoretical predictions. 

In view of this, here we perform quantum dynamical calculation of muon-transfer rates from p\i to C 6+ and 8+ 
and from d/x to 8+ . For this purpose, we use close-coupling approximation to two-component Faddeev-Hahn-type 
dynamical equations jl2|,[l3|]. We are currently investigating muon-transfer rates from muonic hydrogen to other 
heavier elements for a future publication. 
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It is difficult to perform a quantum dynamical calculation of a charge-transfer reaction. In addition, a theoretical 
study of the problems above is extremely complicated due to the large charge of these heavy nuclei and the presence 
of a large number of open channels even at zero incident energy. The large charge of the nuclei leads to a strong 
polarization of the muonic hydrogen in the initial state and a strong final-state Coulomb repulsion. It is difficult to 
incorporate these two effects properly in a dynamical calculation. This is why there are no dynamical calculations of 
these muon-transfer rates. 

Although there is a very large number of open channels in this problem, for a given nuclei the muon is transferred 
predominantly to a few (muonic) atomic labels of the heavy nuclei For example, muon is captured mostly in the 

n = 4 states of C 6+ , and n = 5 states of 8+ . Also these transfers take place mostly to the final muonic-atomic states 
with low angular momenta and transfer rates are negligible for atomic states with angular momenta I > 2. The correct 
dynamical formulation should include all the important transfer channels and we included them in a previous study 
on muon transfer with light nuclear targets [fl4|-jl7|]. However, it is quite impossible now to treat even this reduced 
number of open transfer channels in a dynamicaTcalculation with heavier targets due to convergence difficulties in 
the presence of the large polarization potential and large final-state Coulomb interaction mentioned above. Hence in 
the present treatment we use a two-channel model to calculate transfer to a single final state, where we include the 
elastic and one transfer channels. Different sets of equations are used for the different final states. Eventually, the 
total transfer rate is calculated by summing the different contributions. 

The two-component Faddeev-Hahn-type equation, we use, allows us to introduce explicitly a polarization potential 
in the initial channel and the repulsive Coulomb potential in the final channel. This has the advantage of building 
in the correct asymptotic behavior of the wave function in a low-order close-coupling type approximation ]l8fl . Hence 
as in Ref. [14-hq] we make a two-state close-coupling approximation to the Faddeev-Hahn-type equation in the 
present study ana find that a numerical solution using the present scheme leads to very encouraging agreement with 
experimental transfer rates. 

In Sec. II we present our formalism. Numerical results for muon-transfer rates from muonic hydrogen to carbon 
and oxygen are given in Sec. Ill and compared with those of other investigations. In Sec. IV we present a summary 
and outlook. 



II. THEORETICAL FORMULATION 



The theoretical framework for the present study will be based on the formalism developed in Ref. pq] which was 
used for the study of muon transfer from muonic hydrogen atoms to light charged nuclei, such as, He^~ and Li 3+ . 
Here we shall perform a similar study with heavier charged nuclei, such as, C 6+ and 8+ . The presence of the strong 
Coulomb interaction and the associated large polarization interaction make the present calculational scheme far more 
complicated theoretically and numerically compared to that of Ref. Eq] . The formalism of Ref. |15[ is a generalization 
over that of Ref. |14| for charge-transfer reaction with no final-state Coulomb interaction. In the dynamical equations 
in Ref. [ fl5| the finaTstate Coulomb interaction is explicitly included in the transfer channel. In addition, in the present 
work we explicitly include a polarization potential in the elastic channel. In a coupled-channel approach for atomic 
pro cesses, the coupling to infinite number of p-wave states is responsible for generating the polarization potential 
|L9[ . As it is impossible to include all such states in a numerical scheme, the commonly accepted procedure is to 
replace these coupling terms by a polarization potential as in Ref. This idea has been recently used successfully 
in antiproton- hydrogen and antihydrogen- hydrogen reactions pl| . Next we describe the dynamical equations we use 
based on the close-coupling approximation to Faddeev-Hahn-type two-component equations [ fTi] ]. 

We use units e = h = = 1, where (e) is the muonic mass (charge), and denote, the heavy nuclei (C 6+ or 
8+ ) by 1, the hydrogen isotopes by 2 and muon by 3. Below the three-body breakup threshold, following two-cluster 
asymptotic configurations are possible in the system 123: (23) — 1 and (13) — 2. These two configurations correspond 
to two distinct physical channels, denoted by 1 and 2, respectively. These configurations are determined by the Jacobi 
coordinates (f j3 ,p k ): r 13 = f 3 -rx, p 2 = (r 3 +miri)/(l+mx)-r 2 , r 23 = r 3 -r 2 , pi = (r 3 +m2r2)/(l+m 2 )-ri, 
where ri, m, (i — 1, 2, 3,) are coordinates and masses of the particles respectively. 

Let us introduce the total three-body wave function as a sum of two components 

y(ri,r 2 ,r 3 ) = *i(r 23 ,Pi) + #2(^3, P2), (2.1) 

where ^i(r 23 , p\) is quadratically integrable over the variable f 23 , and ^'2(^13, p 2 ) over the variable f\ 3 . The compo- 
nents ^1 and \l/2 carry the asymptotic boundary condition for channels 1 and 2, respectively. The second component 
is responsible for pure Coulomb interaction in the final state. These components satisfy the following set of two 
coupled equations 

[E -H - V 23 (r 23 ) - C/poiGoi)]*!^, Pi) = [^23(^23) + V 12 {r 12 ) - U c (p 2 )]^ 2 (ri 3 , p 2 ) , (2.2) 
[E-H Q - V 13 (f 13 ) ~ U c {p 2 )]f 2 (r 13 ,p 2 ) = [V 13 (r 13 ) + V l2 {f l2 ) - tfpoi(A)]*i(TWi) > (2-3) 

where E is the center-of-mass energy, Hq is the total kinetic energy operator, Vij(fij) pair potentials (i =/= j = 1, 2, 3), 
Uc is the final-state Coulomb interaction given by 
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Uc{P2 



(Z-1)Z' 
Pi 



(2.4) 



with Z the charge of the heavy nuclei and Z'(= 1) the charge of the hydrogen isotope. Here U po \ is the polarization 
potential given by pf 



Upoi(pi) = --r-T 



92? 
4pf 



Pi > A 



(2.5) 



and z ero other wise. The value of the cut-off parameter A has to be chosen appropriately (see Sec. III). By adding 
Eqs. (2.2) and (2.3) we find that they are equivalent to the Schrodinger equation. 

Distortion potential has been very useful in model and phenomenological description of reaction and scattering in 
nuclear p2| and atomic physics p3| . Although such distortion potentials are unnecessary in a complete solution of the 
Schrodinger equation, they enhance the agreement with experiment in a simplified model description. For example, 
a long-range polarization (distortion) potential has been routinely used in electron-atom scattering |2C],Q. Such a 
distortion potential arising from the polarization o f th e muonic hydrogen isotope due to the bare nuclei is effective in 
the initial channel and has been included in Eq. (^^). This polarization potential for C 6+ or 8+ is much stronger 
and its effect on cross sections much more pronounced than in the case of electron scattering. Hence, for obtaining 
a better agreement with experiment in a model calculation it is prudent to include the polarization potential in the 
elastic channel. There should also be such a polarization potential in the final channel. However, by far the most 
important interaction in the final channel is the Coulomb potential between the proton (or deuteron) and the charged 
muonic atom (A M )( Z ~ 1 ' + . This Coulomb (distortion) potential has also been explicitly included in Eq. ( |2.3| ). This 
will help in obtaining a realistic description of the transfer process as we shall find in the following. 

Because of the strong final-state Coulomb interaction in the present muon-transfer problems it is very difficult to 
develop and s olve succe ssfu lly multichannel models based on Eqs. (^2) and (2^3) above as in Ref. Hence, for 

solving Eqs. ( |2.2[ ) and (2.3) we expand the wave function components in terms of bound states in initial and final 
channels, and project these equations on these bound states. The expansion of the wave function is given by 



*l(^23,Pl) 



fil \pi) R {Z') 



Pi 



(2.6) 



*2(ri 3 ,P2) 



^2^i?if ^(Ifial) {Y £ (p 2 ) ® *Kf 13 )} 00 - 
Pi 



(2.7) 



where nlC are quantum numbers of the three-body final-state, [i\ = 77137712/(7713 + 777 2 ), p 2 = 77737774/ (777,3 + m i)j 
YJm's are the spherical harmonics, R nl L.Qr]) is the radial part of the hydrogen- like bound-state wave function for 

reduced mass pi and charge Z, (pi) and f^x(p2) are ^ ne unknown expansion coefficients. This prescription is 
similar to that adopted in the close-coupling approximation. After a proper angular momentum projection, the set 
of two-coupled integro-differential equations for the unknown expansion functions can be written as 



^ - 2M 1 [/ pol (p 1 ) 



/£Vi)= W(2£ + l) 



/" 

Jo 



dwsinwii^^drasl) 



Z' 

1^23 I 



z 



dpif^Pi) 



nl,fi2 



P1P1C cm Yi m {v 2l TT) I W ^ , 



(2.8) 



dpi pi 

Jo 



z z 



|n3| F121 

-<oo 



Upol{pi 



R[ Z s 'l (\r 23 \)p2PiCZ Yi m (^)/Vi7r . 



(2.9) 



Here k\> = \J2Mi(E - EfJ), k)t' = \J 2M 2 (E - E ( n ') with Mf 1 = m^ 1 + (l + m 2 )- 1 and M 2 X = m^ 1 + {l + m l )- 1 , 

E^ is the binding energy of pair (j3) and gj = AitMj/^ (j = 1, 2), 7 = 1 - mim 2 /((l + tt7i)(1 + 777-2)), ^mL tlic 
Clebsch-Gordon coefficient, u> is the angle between the Jacobi coordinates pi and p 2 , v\ is the angle between 7^3 and 
pi, v 2 is the angle between ri 3 and p 2 . 
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We find that after the projection of the Faddeev-Hahn-type equations (2.2) and (2.3) on the basis states, the initial- 
state pol arization and the final-state Coulomb potentials survive on the left-hand-side of the resultant equations (|2.8| ) 
and (2.9). The presence of the explicit Coulomb potential in the final channel will automatically yield the correct 
physical Coulomb- wave boundary condition in that channel. The explicit inclusion of the polarization potential, 
although has no effect on the boundary con ditio n in t he i nitial channel, substantially improves the the results of the 
truncated model calculation based on Eqs. (2.8) and (2.9) as we shall see in the next section. 

To find unique solution to Eqs. ( |2 . 8[ ) — ( 2 . 9 j , appropriate boundary conditions are to be considered. We impose the 

usual condition of regularity at the origin f[l\o) =0 and /^(O) =0. Also for the present scattering problem with 
1 + (23) as the initial state, in the asymptotic region, two solutions to Eqs. ([2~§1) — (2.9) satisfy the following boundary 
conditions 



P2^ + °0 



Atfcos^pi) , 
A/^^ 1 cos(fcp ) p 2 - J7/2fcf 3 ln2fcf ] p 2 



where K£ l are the appropriate coefficients. For scattering with 2 



(2.10) 

nC/2) , (2.11) 
(13) as the initial state, we have the following 



conditions 



( 2 ) 1„ Ol.( 2 )^ _ 7r £/2) 



fnlM) ~ sin(fcf ) p2-^/2fcf ) ln2^ 2) p 2 

p2 — * + OC 



K% cos(fcf ] p 2 - r ] /2k\ z > \n2k[ z >p 2 - n£/2) , 



(2)- 



.(2). 



(2.12) 
(2.13) 
(2.14) 



where Vi (i — 1, 2) is velocity in channel i. The Coulomb parameter in the second transfer channel is r\ = 2M 2 (Z 



l)/k 



(2) 



The coefficients K£ are obtained from the numerical solution of the Faddeev-Hahn-type equations. The 



cross sections are given by 



4tt 



(K? 2 1 ) 2 



'ls^nl 



fc(l)2 ( D _ 1)2 + (jql + K nly > 

where D = K$K$ - K$K?j. When fc« -► 0: a? s ^ nl ~ l/k[ 1] . The transfer rates are defined by 



(2.15) 



(2.16) 



with vi being the relative velocity of the incident fragments and Nq the liquid-hydrogen density chosen here as 
4.25 x 10 22 cm -3 , because A tr (fc^ 1 - ) — ► 0) ~ const. In our model approach the total muon transfer rate is 



a, 



, Ms 



(2.17) 



III. NUMERICAL RESULTS 



We employ muonic atomic unit: distances are measured in units of a„, where is the radius of muonic hydrogen 
atom. The integro-differential equations are solved by discretizing them into a linear system of equations. The 
) and (|]^) are discretized using the trapezoidal rule and the partial derivatives are discretized 
using a three-point rule |25f| . The discretized equation is subsequently solved by Gauss elimination method. As we 
are concerned with the low-energy limit, only the total angular momentum L — is taken into account. Even at zero 
incident energy, the transfer channels are open and their wave functions are rapidly oscillating Coulomb waves. In 
order to get a converged solution we needed a large number of discretization points. More points are taken near the 
origin where the interaction potentials are large; a smaller number of points are needed at large distances. 

First we solved the system of equations without the polarization potential in the incident channel. However, the 
final-state Coulomb interaction is correctly represented in our model. In this case it is relatively easy to obtain the 
numerical convergence for the system of equations which includes the elastic channel and one transfer channel at a 
time. Finally, the total transfer cross section is calculated by adding the results of different two-channel contributions. 
Without the polarization potential, we needed up to 700 discretization points adequately distributed between and 
50(2^. Near the origin we took up to 60 equally spaced points per unit interval (a p ). 

It was more difficult to obtain convergence with the polarization potential. The polarization potential ( |2.5| ) is taken 
to be zero at small distances below the cut off A. In this case to get numerical convergence we had to integrate to very 



integrals in Eqs. (2.. 
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large distances — up to 300a M . We needed up to 2000 discretization points to obtain convergence. Again we needed 
more points near the origin and less at large distances. For example, near the origin we took up to 60 equally spaced 
points per unit length interval a^; in the intermediate region (p — 10— 20a M ) we took up to 15 equally spaced points 
per unit length interval, and in the asymptotic region (p = 20 — 300a M ) we took up to 5 equally spaced points per unit 
length interval. It is well-known that the results for the cross section is sensitive to the value of the cut off A of the 
polarization potential. The short-range potential of the present problem extends to about 25a M . We considered the 
polarization potential in the asymptotic region p\ > A ~ 75a^. For a small variation of A in this region from 75a^ to 
about 120(2^, we find the transfer cross sections to be reasonably constant and the reported transfer cross sections of 
this study are the averages of these cross sections. If A is increased past 120(2^, the effect of the polarization potential 
on the cross sections gradually decreases and finally disappears. If A is decreased much below 75a M , the cross sections 
become rapidly varying function of A and could become unphysically large. The range of A values (here between 75a^ 
and 120a M ) for which the cross sections are slowly varying smooth functions should increase with the charge of the 
bare nucleus. 

We present partial muon-transfer rates A*^ and total transfer rates A^ t calculated using the formulation of last 
section. In this work using the model of Sec. II we calculate the low-energy muon-transfer rates from (pp)i s to C 6+ 
and 8+ and from (dp)i s to 8+ . From other theoretical |^,[llj investigations it was concluded that in the case of 
C 6+ the transfer takes place predominantly to the n = 4 state and for 8+ it happens to the n — 5 state, which is 
also found to be true in our model calculation. Hence in this work we only present rates for the I = 0, 1, 2 states of 
the n — 4 and 5 orbitals of carbon and oxygen, respectively. The contribution of the higher angular momentum states 
to the total transfer cross section is very small. Numerically converged results were obtained in these cases. The low 
energy partial rates A^_ >ni /10 10 sec -1 and total rates A*J, t /10 10 sec" 1 are presented in Tables I, II, and III together 
with the results of other theoretical and experimental works. 

First we comment on the results in Table I for muon transfer from (pp)i s to C 6+ . The partial transfer rates without 
the polarization potential increases with decreasing center-of-mass energy E and saturates to a constant value for 
E < 0.04 eV for 4s, 4p, and 4d states of muonic carbon. In the case of C 6+ the transfer takes place predominantly to 
the 4s state. These qualitative behaviors are also true after the inclusion of the polarization potential and were also 
true in a previous theoretical study. However, after the inclusion of the polarization potential the 4s and 4p transfer 
rates are enhanced by about a factor of 1.5. The present total transfer rate of 8.5 x 10 10 sec -1 is about three times 
larger than the previous theoretical calculation of 2.8 x 10 10 sec -1 |ll|. We quote two experimental results in this 
case: (9.5 ± 0.5) x 10 10 sec -1 Q and (5.1 ± 1.0) x 10 10 sec -1 [Q. The present theoretical result lies in between these 
two somewhat conflicting experimental results. 

In the case of muon transfer from (pp)i s to 8+ , we find from Table II that transfer takes place predominantly to 
the 5s and 5p states of muonic oxygen. Again the transfer rates saturate and attain constant values for E < 0.04 eV. 
The transfer rate is higher in the 5s state and lowest in the 5d state. This behavior remain true after the inclusion of 
the polarization potential, when the transfer rate to the 5p state increases by a factor of more than two whereas the 
contribution to the 5s state increases by a factor of 1.5. The present total transfer rate of (7.7 ± 0.5) x 10 10 sec -1 is 
about 1.5 times larger than the previous theoretical calculation of 5.6 x 10 10 sec -1 |11| and in reasonable agreement 
with the recent experimental rate of (8.5 ± 0.2) x 10 10 sec" 1 [Q. 

Finally, in the case of muon transfer from (dp)i s to 8+ , we find from Table III that transfer also takes place 
predominantly to the 5s state of muonic oxygen. The contribution to the transfer rate due to 5s state is two times 
as large as the contribution due to the 5p state. Again the transfer rate saturate and attain a constant value for 
E < 0.04 eV. After the inclusion of the polarization potential the transfer rate to the 5s state increases by a factor of 
1.5. The present total transfer rate of (4.4 ± 0.6) x 10 10 sec -1 is in reasonable agreement with the recent experimental 
rate of 5.5 x 10 10 sec" 1 §. 

The Z dependence of the transfer rates from a specific hydrogen isotope to X z+ has been a subject of interest. 
Although for large Z these rates increase linearly with Z, there is no general behavior for small Z. The most 
recent experimental transfer rates decrease when we move from the system pp — C 6+ to pp — 8+ MM- Through 
our dynamical calculation we have been able to reproduce this behavior. Our calculation is also consistent with the 
experimentally observed isotope effect, e.g., the transfer rate decreases when we move from pp — 8+ to dp — 8+ 



IV. CONCLUSION 



We have studied muon-transfer reactions from muonic hydrogen to carbon and oxygen nuclei employing a full 
quantum-mechanical few-body description of rearrangement scattering by solving the Faddeev-Hahn-type equations 
using close-coupling approximation. To provide the correct asymptotic form of the wave function in the transfer 
channel, the final-state Coulomb interaction has been incorporated directly into the equations. We also included a 
polarization potential at large distances in the initial channel. It is shown that in the present approach, the application 
of a close-coupling-type ansatz leads to satisfactory results for direct muon-transfer reactions from muonic hydrogen 
to C 6+ and 8+ . In the case of muon transfer from (pp)i s to C 6+ the present transfer rate to the 4s state of muonic 
carbon is about 1.5 times larger than that to the 4p state. For muon transfer from (pp)i s to 8+ , the present 
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transfer rate to the 5s state of oxygen is about twice as large as that to the 5p state. The present rates are much 
larger by factors of about two to three compared to the calculation of Ref . . Finally, in the case of muon transfer 
from (<i/i)i s to 8+ , the present transfer rate to the 5s state is large compared to that to the 5p state. In all cases 
the inclusion of the polarization potential improves the agreement with experiment and our final transfer rates are 
in encouraging agreement with recent experiments 0-0]. The present rates for oxygen from (p/J.)i s and (d/i)i s are 
in agreement with the observed isotope effect ||: the transfer rate increases with the decrease of the mass of the 
hydrogen isotope. Because of the present promising results for the muon-transfer rates for Z = 6 and Z = 8 it seems 
useful to make future applications of the present formulation for larger targets. Calculations involving nuclei of higher 
charges (Ne 10+ , S 16+ , Ar 18+ etc.) are in progress. 
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Table I. Low energy partial A' r s ^„;/10 10 sec _1 and total A*o t /10 10 sec _1 muon transfer rates 
reduced to liquid-hydrogen density Ao = 4.25 x 10 22 cm~ 3 from muonic hydrogen (pfi)i s to 
hydrogen-like excited state of muonic carbon (C M )j£t 4 . 
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Table II. Low energy partial A5 1 s _,„;/10 10 sec _1 and total A*o t /10 10 sec _1 muon transfer rates 
reduced to liquid-hydrogen density iVo = 4.25 x 10 22 cm -3 from muonic hydrogen (p/^)i s to 
hydrogen-like excited state of muonic oxygen (O ll ) r ^L 5 . 
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Table III. Low energy partial A' r s _ +ni /10 10 sec _1 and total A*Q t /10 10 sec" 1 muon transfer rates 
reduced to liquid-hydrogen density Nq — 4.25 x 10 22 cm~ 3 from muonic hydrogen (dfi)i s to 
hydrogen-like excited state of muonic oxygen (0 /J )^ 5 . 
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